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Abstract 

Examples of test electromagnetic waves on a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) background are constructed from explicit 
perturbations of the FLRW space-times describing gravitational waves 
propagating in the isotropic universes. A possible physical mechanism 
for the production of the test electromagnetic waves is shown to be 
the coupling of the gravitational waves with a test magnetic field, 
confirming the observation of Marklund, Dunsby and Brodin [Phys. 
Rev. D62, 101501(R) (2000)]. 
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1 Introduction 



Marklund, Dunsby and Brodin [lj have made the important observation that 
gravitational wave perturbations of Friedmann-Lemaitre-Robertson- Walker 
(FLRW) cosmological models coupled to weak magnetic test fields can gener- 
ate electromagnetic waves. The magnitude of this effect has been estimated 
in pp. They have demonstrated this phenomenon using the gauge-invariant 
and covariant perturbation theory of Ellis and Bruni [2]. This theory has also 
been used to construct explicit solutions of the Ellis-Bruni perturbation equa- 
tions describing gravitational waves propagating in FLRW universes [3] (see 
also [1]) from a point of view which differs significantly from that of [lj. The 
histories of the wave fronts of the latter waves are particularly simple fami- 
lies of null hypersurfaces which arise naturally in the FLRW space-times. In 
this paper we demonstrate how test electromagnetic waves can be constructed 
from these explicit gravitational waves and how these electromagnetic waves 
can be viewed as arising from the interaction of the gravitational waves with a 
test magnetic field, thereby supporting the observation of Marklund, Dunsby 
and Brodin. 

We utilize a test magnetic field which, by its nature, does not perturb the 
isotropic cosmology but is nevertheless a violation of isotropy. The gravita- 
tional waves we use do perturb the isotropic cosmology and this perturbation 
also breaks the isotropic symmetry since the gravitational waves are uni- 
directional and thus at each point of the isotropic cosmological model their 
histories have a unique null propagation direction. The behavior of magnetic 
fields in cosmological models has been extensively and carefully studied from 
diverse physical viewpoints in [5] — |10j . 

The gravitational waves which we utilize in this paper correspond to the 
simplest (from a geometrical point of view) type of gravitational radiation 
which can propagate in isotropic cosmologies. The histories of their wave 
fronts are naturally occurring null hypersurfaces in the isotropic cosmological 
models and their propagation direction in these space-times is null, geodesic 
and shear-free. It is thus of some interest to examine the Marklund, Dunsby 
and Brodin observation in terms of them. As a measure of the strength 
of the gravitational waves used in this paper, and of the electromagnetic 
waves which result from their interaction with a test magnetic field, we find 
that the gravitational field of the gravitational waves (the perturbed Weyl 
tensor) is proportional to R~ 2 (see eq.(2.34)) where R(t) is the scale factor 
of the isotropic universe, while the electromagnetic field (see eq. (3.17)) is 
also proportional to R~ 2 . The coefficients of R~ 2 in both cases involve an 
arbitrary analytic function whose appearance is a characteristic of electro- 
magnetic radiation (which is shear-free [H] in the optical sense that the 
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null propagation direction in space-time, given by the gradient of the func- 
tion <p in (2.28) below, is shear-free) and of shear-free (in the optical sense) 
gravitational radiation [12J. 

The interaction of gravitational waves and electromagnetic fields has led 
to proposals for a mechanism to detect gravitational waves [13] [15]. Our 
approach is influenced by the fundamental paper by Szekeres [TU] on the in- 
teraction of gravitational waves with matter and with electromagnetic waves. 
With the inclusion of matter many modelling possibilities open up. For ex- 
ample, a recent model for the generation of gravitational waves from matter 
and electromagnetic waves can be found in [TT] . 

This paper is organized in the following way: in section 2 the relevant 
Ellis-Bruni perturbation equations (for tensor perturbations) are listed and 
the explicit solutions derived in [3] are summarized. Analogous test electro- 
magnetic waves having the same wave fronts as these gravitational waves will 
be required and they are described in section 3. In section 4 the construction 
of test electromagnetic waves of the type considered in section 3 are derived 
from the gravitational waves of section 2. In addition in this section these 
electromagnetic waves are shown to be capable of being interpreted as re- 
sulting from the interaction of the gravitational waves with an explicit test 
magnetic field. The paper ends with a brief discussion in section 5. 

2 Gravity Waves in FLRW Universes 

We make use of a four dimensional space-time manifold with a metric tensor 
having components g ab in a local coordinate system {x a }. This manifold con- 
tains a preferred congruence of time-like world-lines which are the integral 
curves of a vector field having components u a in the coordinate system {x a } 
and which satisfies u a u a = —1. The electric and magnetic parts of the Weyl 
tensor C a bcd are defined respectively by 

E ab = C apbq u p u q , H ab = C: pbq u p u q . (2.1) 

These are equivalent to C a bcd and the star superscript denotes the dual 
C a P bq = lVa P rs C rsbq with r] abcd = y/=g e abcd , g = det(g ab ) and e abcd the Levi- 
Civita permutation symbol. We note that for the Weyl tensor the left and 
right duals are equal. The energy-momentum-stress tensor T ab = T ba de- 
scribing the matter distribution can be decomposed with the respect to the 
vector field u a to read [18J 

T ab = fiu a u b + ph ab + q a u b + q b u a + n ab , (2.2) 
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with fx the density of matter measured by the observer with 4- velocity u a , 
h ab = g ab + u a u b the projection tensor, p the isotropic pressure, q a the en- 
ergy flow measured by the observer with 4-velocity u a and 7r ab = n ba the 
anisotropic stress. Here 

q a u a = 0, n ab u b = 0, 7r a a = 0. (2.3) 

The covariant derivative u a - b is decomposed into 

1 

u a -b = u ab + a ab + -9 h ab -u a u b , (2.4) 

with the dot in general indicating covariant differentiation in the direction of 
u a (and thus in particular ii a = u a - b u b ). Here 

Uab = W[o;6] + U [a U b \ , (2.5) 

is the vorticity tensor (with u ab = —uj ba , u ab u b = and square brackets 
denote skew symmetrization) , 

Oab = U(a-b) + U{a Ub) - 7^ Kb , (2.6) 

is the shear tensor (with a ab = a ba , a a a = , a ab u b = and round brackets 
denote symmetrization) and 

= u a , a , (2.7) 

is the expansion or contraction scalar. 

In the isotropic FLRW space-times with u a the 4-velocity of matter we 
have q a = and n ab = and thus ( 12.2ft specialises to the perfect fluid form 

T ab = fiu a u b + ph ab , (2.8) 

with, in addition, h b p^ b = and h b a fi tb = 0. We also have in this case u a = 0, 
u ab = 0, a ab = and h b 9 jb = which has the effect of simplifying (12 .4p to 

u a - b =\0h ab . (2.9) 
In this case 9 satisfies the simplified Raychaudhuri equation 

o + \e 2 = -\^ + zp) , (2.io) 

and the equations T ab . b = reduce in this case to the single equation 

fi + 9(fi + p) =0 . (2.11) 
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The space-times are now necessarily conformally flat and thus E ab = 
and H ab = 0. The metric tensor g ab is given via the line-element in the 
Robertson- Walker form 

d s^RW {dXl)2 + {dx2)2 X {dx3)2] -dt\ (2.12) 
1 + |r 2NZ 

with R(t) the scale factor, r 2 = (x 1 ) 2 + (a; 2 ) 2 + (x 3 ) 2 and k(= 0, ±1) is 
the Gaussian curvature of the t = constant space-like hypersurfaces. In 
these coordinates u a d/dx a = d/dt and x a (a = 1,2,3) are constant on each 
integral curve of the vector field d/dt. In addition 6 = 9(t), \i = fi(t) and 
p = p(t) are obtained from the scale factor R{t) as 6 = 3R/R with fi, p 
given by ( 12. lip along with an equation of state while R(t) is derived from 
the Friedmann equation (whose time derivative in this case coincides with 
(12111 : see, for example [T8]). 

In the Ellis-Bruni [2] theory perturbations of the FLRW models in partic- 
ular are described by the basic gauge invariant variables E ab , H ab , a ab , u ab , it a , 
h b a tb , h b a fi tb , h b a p tb , q a and n ab , which are taken to be small of first order. 
The differential equations determining these variables are obtained by retain- 
ing first order small terms in the Ricci identities, the Bianchi identities and 
the equations T ab . b = with T ab given by (12.21) . There are two exceptions to 
this procedure: the propagation equation for 6 along the integral curves of u a 
(Raychaudhuri's equation), which follows from the Ricci identities, and the 
propagation equation for /i along the integral curves of u a , which follows from 
T ab . b = 0. These equations are converted into differential equations for the 
gauge invariant variables listed above by projecting their gradients orthogo- 
nal to u a (see [2]). For perturbations that exclusively describe gravitational 
waves propagating through FLRW universes we only require the variables 
E ab , H ab , a ab , n ab (the so-called "tensor" perturbations), with the remain- 
ing gauge invariant variables vanishing. The linear equations determining 
these first order quantities are [3] 

2 1 

<?ab + g a ab - -7T ab + E ab = , (2.13) 

and 

H ab = ~(T (a 9]C r] b)fgc u f , (2.14) 
from the Ricci identities and 

E bd -, d = -l^ d ;d, (2.15) 
H bd ;d = 0, (2.16) 
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E bt + 6E bt = -u r Hf, r]^ rsd --n bt --6TT b \ (2.17) 

2 6 

H bt + 9H bt = u r E^ d rfy sd -^ b rad ^ d u r , (2.18) 

from the Bianchi identities. To obtain simple solutions of these equations we 
look for solutions which have an arbitrary dependence on a scalar function 
(j){x a ) by writing 

°ab = Sab F{<f>) , 7l ab = U ab F{<f)) . (2. 19) 

We will then substitute these into (12.131) and ( 12. 14ft to obtain expressions 
for E a b and H ab which depend on F and its derivative with respect to 4>. 
These, along with (12.191) . are then substituted into (2.15)-(2.18) to arrive at 
differential equations for s a b and U ab . We note that 

Sab = S ba , S ab U b = 0, g* S ab = s\ = 0, (2.20) 

with similar equations holding for U a b- When the substitutions indicated 
above are carried out the consistency of the resulting equations requires [3] 

g ab ^,a <P,b = 0, (2.21) 

(the comma denoting partial differentiation with respect to x a ) so that the 
hypersurfaces <fi(x a ) = constant are null, and the following consistent equa- 
tions (also consistent with (I2.20p ) must hold: 

s ab 0, b = O, rP fe 0, fe = O, (2.22) 

along with 

s ab , b = , U ab . b = , (2.23) 

the propagation equation for s ab along the generators of the null hypersurfaces 
4>(x a ) = constant, 

s ab '' c <P,c + fy 4 ; d - \e^j s ab = - l -j>W b , (2.24) 
and the wave equation, 

ab;d _ l e -ab + (_q2 + 3 \ g ab = _ jjoft _ 2 _Q ^ab ^ ^ 

'3 V 3 T 6V 3 y J 

The null hypersurfaces <p(x a ) = constant are the histories in the FLRW 
space-times of the wave fronts of the gravitational waves described by these 
perturbations (see [3], [1]). To exhibit explicit solutions we first choose these 
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null hypersurfaces. To obtain some naturally occuring null hypersurfaces 
which will lead to surveyable solutions we start by writing the Robertson- 
Walker line-element (12.121) in the form [19] 

ds 2 = R 2 (t){dx 2 + Po 2 f 2 (dy 2 + dz 2 )} - dt 2 , (2.26) 

with p = 1 + (K/4)(y 2 + z 2 ), K = constant, / = f(x). The following cases 
are allowed: (i) if k = +1 then K = +1 and f(x) = sinx; (ii) if k = then 
K = 0, +1 with f(x) = 1 when K = and f(x) = x when K = +1; (hi) if 
k = —1 then K = 0, ±1 with f(x) = ~ e x when K = 0, f(x) = sinhx when 
K = +1 and /(x) = coshx when K = — 1. The details of these special cases 
are given in [3]. The following equations are satisfied in the cases (i)-(iii): 

f" + kf = 0, {ff + kf 2 = K 1 (2.27) 

with the prime denoting differentiation with respect to x. A convenient 
family of null hypersurfaces is given by 

(j)(x a ) =x- T{t) = constant , (2.28) 

with dT/dt = R- 1 . 

Using the null hypersurfaces (12.281) and the equations (2.23)-(2.25) for 

s ab , U ab subject to (2.20) (and similar conditions on IT 1 '') and (2.22) yields 
solutions in the form 

s ab = sm a m b + sm a rh b , (2.29) 

U ab = n m a m 6 + n?n a m fr , (2.30) 

with m a given by the 1-form m a dx a = Rp^ 1 f dC,/ \/2 with ( = y + iz (and 
thus m a m a = = m a rh a and m a m a = +1) and the bar denoting complex 
conjugation. In addition 

,l "g({,x,t), (2.31) 



Rf 



fl = -^&(Dg + RQ) . (2.32) 

R 2 j 

Here Q is a complex analytic function of its argument and D = d/dx+R 8/ dt. 
Also Q satisfies 

D 2 g + kg = o, (2.33) 

with k = 0, ±1 and this can easily be solved with each case involving two 
arbitrary complex analytic functions, of ( and x — T = 0, of integration. The 
corresponding perturbed Weyl tensor is given by 

E ab + lR ab = _M JL(£ F (0)) m a m * j ( 2 .34) 
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which is a Petrov type N Weyl tensor with degenerate principal null direction 
<f) t a, confirming the interpretation of the perturbations of the FLRW space- 
times described here as being due to gravitational waves having the null 
hypersurfaces <fi = constant as the histories of their wave fronts. Further 
properties of these waves are discussed in [3]. 

3 Electromagnetic Waves in FLRW Universes 

Electromagnetic test fields on the space-time described at the beginning of 
section 2 are encoded in a skew-symmetric tensor with components F ab = 
—F ba having electric and magnetic parts defined by [IB] 

E a = F ab u b , H a = F: b u b , (3.1) 

with F* b = \rj ab rs F rs . The vectors E a , H a are equivalent to a knowledge of 
F ab with 

F ab = u a E b - u b E a - r] abcd u c H d . (3.2) 

If these are source-free electromagnetic test fields on the FLRW space-times 
then they must satisfy Maxwell's equations in the form [18] 

E a , a = , H a ;a = , (3.3) 

and 

E a + ^6E a = - V abed u b H e , d , (3.4) 

H a + hH a = V abed u b E e . d . (3.5) 

To obtain solutions analogous to the gravitational waves described in section 
2 we first introduce a 4-potential o a with [20] 

a a u a = 0, a a . a = 0, (3.6) 

and 

Kb = Vb;a ~ & a;b ■ (3.7) 

Using the Ricci identities along with the conformal flatness of the FLRW 
space-times we have 

2 

o"a;dc - o a . cd = -fi g a[c a d] + {fi + p) u a U[ c a d] . (3.8) 
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This helps to establish that (13.31) are now satisfied automatically. In addition 
( 13.41) reduces to the wave equation 

v a ' 4 ,d = \(fi-p)v a , (3.9) 

and (3.5) is automatically satisfied. The equations to be satisfied by a a are 
(13.61) and (13. 9p . Following (12.191) we look for solutions of the form 

o a = s a F(<f>) , (3.10) 

with F an arbitrary function of <fi(x a ). With s" ^ we arrive again at (I2.2ip 
along with 

s a , a = } s a <f), a = 0. (3.11) 

We also obtain the propagation equation for s a along the generators of the 
null hypersurfaces (p(x a ) = constant, 

s a -> b ^ b + ^/s a = 0, (3.12) 

and the wave equation, 

s a - 4 . 4 =\^-p)s a . (3.13) 

Solving these with <p(x a ) given by ( 12.281) we find 

s a = s m a + s rh a , (3.14) 
with m a given following ( I2.29P and 

so = ^jF{C,,x,t) , (3.15) 

with T a complex analytic function required to satisfy 

DF = 0, (3.16) 

with the operator D defined following (2.32). Thus T = J-'iC, x — T) and is 
otherwise arbitrary. The corresponding electric and magnetic fields (13.11) are 
given via 

E a +iH a = ^- ^{FF) m a , (3.17) 

analogous to (I2.34p . The radiative nature of this electromagnetic field, with 
propagation direction is evident from ( 13.171) . 
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4 Electromagnetic Waves from Gravity Waves 



We now demonstrate how electromagnetic waves of the type described in 
section 3 can be constructed from the gravitational waves given in section 
2. The construction is not unique (see the discussion in the next section) 
but it can be given the interpretation of arising from the coupling of the 
gravitational waves of section 2 with a test magnetic field in agreement with 

PQ- 

The gravitational waves of section 2 are obtained from a complex analytic 
function Q((, x, t) which satisfies (I2.33p . We see from the argument in section 
3 that the electromagnetic waves given there are obtained from a complex 
analytic function JF(£,x, t) which satisfies (13. 16j) . Given Q satisfying (I2.33|) 
and f{x) satisfying (12.271) we define 

r{C,x,t) = fDg-fg. (4.i) 

Clearly this function satisfies (13.161) and so T = ^(C, x—T). With this choice 
of T electromagnetic waves in section 3 corresponding to the gravitational 
waves of section 2 are described by the 4-potential (13.101) with s a given by 
dSH and fl3~m 

A source-free test magnetic field on the FLRW space-times must satisfy 
(3.3)-(3.5) with H a ^ and E a — 0. A simple example of such a field is 
given by 

#a=^, K = hU, b - (4.2) 

This field can be derived from a 4-potential in the manner of (13. 7\} with (13.61) . 
In this case the potential is given by the 1-form 

, „ y dz — z dy , 
H a a dx a = y — y - , 4.3 

or equivalently by the covariant vector 

i((m n -(m a ) 
H(T = 2V2Rf ' (4 ' 4) 

with m a given following (2.30). We note in passing that any source-free 
magnetic test field on the FLRW space-times is given by H a = R~ 1 q ja with 
q a u a = and g ab q a b = and (14.21) corresponds to q = q(x) given by dq/dx = 

Observing that the function JF in (14.11) is linear in the arbitrary analytic 
function Q and its derivative DQ we first note from (2.29)-(2.32) that 

2 2 

s ab = -^-Gm a m b -^-Gm a m b , (4.5) 
Rj Rf 
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and 



U ab --6s ab = -^-Dgm a m b -^-Dgm a m b . (4.6) 



3 R 2 f R 2 f 

Hence using ( 14.4ft we arrive at 



- 2 p^R f f s ab H a b - Po l R 2 f 2 (u ab - 2 -6 s ab ^j H a b = (Jm a + J rh a ) , 

(4.7) 

with 

J = ^, («) 

and T is given by ( 14. ip . The right hand side of ( 14.71) is the 4-potential of a 
radiative Maxwell test field on the FLRW space-times of the type described 
in section 3. The left hand side of ( 14. Tl) provides a Marklund-Dunsby-Brodin 
physical interpretation of the origin of this electromagnetic field as a coupling 
of the gravitational waves described by s ab and U ab with the magnetic test 
field described by the 4-potential #cr a , due to the appearance of the products 
s ab H a b and U ab H a b . 



5 Discussion 

We make two observations regarding the electromagnetic waves constructed 
in section 4 from the gravitational waves in section 2. The first observation 
concerns the lack of uniqueness of the resulting electromagnetic waves. To 
see this we note that if Q and TC are two complex analytic functions satisfying 
(JOSDthen 

F((,x,t) = kHg + DHDG , (5.1) 

satisfies ( 13.161) . In view of ( 12.271) a choice of H is simply H = f(x). In this 
case 

F(C,x,t) = kfg + fDg. (5.2) 

It is easy to see from (13. 1 7[) that these waves are distinct from those ob- 
tained using T given in (14. ip . In addition we note that there are no waves 
of this type for the case (ii) with k = K = given following (12. 26j) . in con- 
tradistinction to the example ( 14. ip . On the other hand since ( 15.11) is a linear 
combination of g and Dg these electromagnetic waves also lend themselves 
to the interpretation of arising from the coupling of the gravitational waves 
of section 2 to the test magnetic field given in section 4. 

The significance of utilizing a test magnetic field rather than a test electric 
field in (14. 21) is clarified by the following observation: A source-free test 
electric field given by (3.1), (3.6) and (3.7) with H a = has a potential 
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1-form of the form E a a dx a = a± dx + a 2 dy + a 3 dz with a a (a = 1,2,3) 
functions of x,y,z,t. Maxwell's equations (3.3)-(3.5) with H a = imply 
that E& a = -R _1 /(x, y, z) u a , for some function /, up to a gauge transformation. 
The scalar products of this 4-potential with s ab and U ab are therefore zero. 

Finally we remark that astrophysical gravitational waves are generally 
low frequency and in universes containing ionized matter the correspond- 
ing electromagnetic waves propagate poorly. In our treatment no mention 
has been made of the frequency of the resulting electromagnetic waves with 
4-potential (4.7). This would be introduced via a Fourier analysis of the ar- 
bitrary analytic function Q. To further study this electromagnetic radiation 
and follow the dissipative consequences of its interaction with the matter 
and gravitational radiation would require calculating its perturbative effect 
on the cosmological model. In this paper this electromagnetic radiation is 
weak but is still considered a test field. 
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